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Gauge Invariance and Confinement in Noncompact Simulations of SU(2) 

Kevin CahiU'fl 

'^Department of Physics and Astronomy, University of New Mexico, 
Albuquerque, New Mexico 87131-1156, U. S. A. 

Wilson loops have been measured at strong coupling, f3 = 0.5, on a 12* lattice in a noncompact simulation 
of pure SU(2) in which random compact gauge transformations impose a kind of lattice gauge invariance. The 
Wilson loops suggest a confining potential. 



1. INTRODUCTION 

In 1980 Creutz ||l[ displayed quark confine- 
ment at moderate coupling in lattice simulations 
of both abelian and nonabelian gauge theories. 
Whether nonabelian confinement is as much an 
artifact of Wilson's action as is abelian confine- 
ment remains unclear. The basic variables of Wil- 
son's formulation are elements of a compact 
group and enter the action only through traces 
of their products. The Wilson action has false 
vacua H] which affect the string tension |^J^ . 

In simulations of SU(2) with gauge-invariant 
potential barriers between the true vacuum and 
the false vacua, the string tension has been seen 
to drop Q or even vanish Q. 

To examine these questions, some physicists 
have introduced lattice actions that are noncom- 
pact discretizations of the continuum action with 
fields as the basic variables. For U(l) these non- 
compact formulations are accurate for all cou- 
pling strengths for SU(2) they agree well with 
perturbation theory at very weak coupling Q . 

This report relates the results of measuring 
Wilson loops at strong couphng, f3 = A/g"^ = 
0.5, on a 12'^ lattice in a noncompact simulation 
of SU(2) gauge theory without gauge fixing or 
fermions. In this simulation the fields are sub- 
jected to random compact gauge transformations 
which restore a semblance of lattice gauge invari- 
ance. 



2. NONCOMPACT METHODS 

Patrascioiu et al. performed the first noncom- 
pact simulations of SU(2) by discretizing the clas- 
sical action and fixing the gauge They saw a 
Coulomb force. 

Later simulations |^ were carried out with an 
action free of spurious zero modes, for which it 
was not necessary to fix the gauge. The Wilson 
loops of these simulations showed no sign of quark 
confinement. A possible explanation of this neg- 
ative result is that noncompact actions lack an 
exact lattice gauge invariance. Yet if one sub- 
jects the fields to random compact gauge trans- 
formations during each sweep, then one may be 
able to restore a kind of gauge invariance to the 
simulation [||. Here I report the results of such 
a gauge-invariant simulation in which the Wilson 
loops fall off exponentially with the area of the 
loop. 

In both the earlier simulations without gauge 
invariance and the new simulation with gauge in- 
variance, the fields are constant on the links of 
length a, the lattice spacing, but are interpolated 
linearly throughout the plaquettes. In the plaque- 
tte with vertices n, n + e^, n + e^, and n + e^ + e^, 
the field is 



(- 



?i^)A°(n + e„) 



+ {n^ + l-^)Al{n), 
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and the field strength is 



(1) 



(2) 



2 



The action S is the sum over ah plaquettes of the 
integral over each plaquette of the squared field 
strength, 



(3) 



The mean value in the vacuum of a euclidean- 
time-ordered operator Q is approximated by a 
ratio of multiple integrals over the Aj^(n)'s 



(TQ(A))o 



-S(A) 



(4) 



which one may compute numerically. Macsyma 
was used to write most of the Fortran code 
for the present simulation. 

3. GAUGE INVARIANCE 

To restore gauge invariance, the fields are sub- 
jected to random compact gauge transformations 
during every sweep, except those devoted exclu- 
sively to measurements. At each vertex n a ran- 
dom number r is generated uniformly on the in- 
terval (0, 1); and if r is less than a fixed proba- 
bility, set equal to 0.5 in this work, then a ran- 
dom group element U{rL) is picked from the group 
SU{2). The fields on the four links coming out of 
the vertex n are then subjected to the compact 
gauge transformation 



(5) 



and those on the links entering the vertex to the 
transformation 



-igaA''^(n-e^)Tb 



U{n)e 



-igaA'Un-ef^)Ta 



(6) 



4. WILSON LOOPS 



The quantity normally used to study confine- 
ment in quarkless gauge theories is the Wilson 
loop W(r, t) which is the mean value in the vac- 
uum of the trace of a path-and-time-ordered ex- 
ponential of a line integral of the connection 
around an r x i rectangle 



W{r,t) = [l/d) (tr7'Te~'^i''^^^"''"-)o 



where d is the dimension of the generators T^. 
Although Wilson loops vanish in the exact the- 
ory , Creutz ratios xl'', i) of Wilson loops de- 
fined [|2) as double differences of logarithms of 
Wilson loops are finite. For large t, xi'f'^t) ap- 
proximates (a^ times) the force between a quark 
and an antiquark separated by the distance r. 

In this simulation the data are not yet suffi- 
cient to allow one to determine the Creutz ratios 
beyond the 3x4 loop. The Wilson loops there- 
fore have been fitted to an expression involving 
Coulomb, perimeter, scale, and area terms. 

5. MEASUREMENTS AND RESULTS 

It will be useful to compare this simulation with 
an earlier one |^ in which the fields were not sub- 
jected to random gauge transformations. Both 
simulations were done on a 12^ periodic lattice 
with a heat bath. The earlier simulation consisted 
of 20 independent runs with cold starts. The first 
run had 25,000 thermalizing sweeps at inverse 
couphng /? = 2 followed by 5000 at /? = 0.5; the 
other nineteen runs began at /3 = 0.5 with 20,000 
thermalizing sweeps. There were 59,640 Parisi- 
assisted measurements, 20 sweeps apart. 

Noncompact Wilson loops at (3 ~ 0.5 



r X i 

a a 


Not invariant 


Invariant 


1 x 1 


0.402330(6) 


0.254564(8) 


2x2 


0.085426(4) 


0.018711(6) 


3x3 


0.018080(2) 


0.001429(3) 


4x4 


0.003993(1) 


0.000117(3) 


5x5 


0.000893(1) 


0.000014(6) 


6x6 


0.000201(0) 


0.000004(2) 



(7) 



The present simulation with random gauge 
transformations is very noisy. So far it consists of 
13 runs, all with cold starts and 20,000 thermal- 
izing sweeps. Wilson loops have been measured 
every five sweeps for a total of 689,684 measure- 
ments. The values of the diagonal Wilson loops so 
obtained are listed in the table. The errors have 
been estimated by the jackknife method, with all 
measurements in bins of 100 considered to be in- 
dependent. 
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Figure 1. The negative logarithms of Wilson 
loops with the perimeter factor canceled are plot- 
ted against the area rt of the loop. The loops of 
the gauge-invariant simulation are represented by 
bullets; those of the earlier simulation by circles. 



The Wilson loops of the gauge-invariant sim- 
ulation fall off much faster with increasing loop 
size than do those of the earlier simulation. Be- 
cause the data do not accurately determine all 

the Creutz ratios, I have fitted both sets of loops, 
including the non-diagonal loops, to the formula 



W{r, t) 



^a+b(r/t+t/r)-2c {r-\-t)-drt 



(8) 



in which a is a scale factor, h a Coulomb term, c a 
perimeter term, and d an area term. For the sim- 
ulation without random gauge transformations, I 
found a « 0.25, h w 0.20, c « 0.39, and d w 0.00. 
For the simulation with random gauge transfor- 
mations, I found a ~ 0.60, h » 0.20, c » 0.56, 
and d « 0.11. In the gauge- invariant simulation, 
the coefficient of the area-law term is about two 
orders of magnitude larger than in the earlier sim- 
ulation which lacked gauge invariancc. 

To exhibit the rcnormalizcd quark-antiquark 
potential, I have plotted in the figure the negative 
logarithms — \ogi^ (e^^*^''+*'T4^(r, t)) of the Wilson 
loops with the perimeter term removed. Apart 
from the uncertain value of W^(6,6), the bigger 
loops of the gauge-invariant simulation, repre- 



sented by bullets, display an area law; whereas 
the larger loops of the earlier simulation, repre- 
sented by circles, show an essentially flat poten- 
tial. The smallest loops reflect the symmetrized 
Coulomb term oc {t/r + r/t). 
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